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ABSTRACT 

Future weak lensing measurements of cosmic shear will reach such high accuracy that second order 
effects in weak lensing modeling, like the influence of baryons on structure formation, become impor¬ 
tant. We use a controlled set of high-resolution cosmological simulations to quantify this effect by 
comparing pure N-body dark matter runs with corresponding hydrodynamical simulations, carried 
out both in non-radiative, and in dissipative form with cooling and star formation. In both hydro- 
dynamical simulations, the clustering of the gas is suppressed while that of dark matter is boosted 
at scales fc > 1 hMpc“^. Despite this counterbalance between dark matter and gas, the clustering of 
the total matter is suppressed by up to 1 percent at 1 < fc < lO/iMpc”^, while for k ss 20/iMpc~^ it 
is boosted, up to 2 percent in the non-radiative run and 10 percent in the run with star formation. 

The stellar mass formed in the latter is highly biased relative to the dark matter in the pure N-body 
simulation. Using our power spectrum measurements to predict the effect of baryons on the weak 
lensing signal at 100 < I < 10000, we find that baryons may change the lensing power spectrum by 
less than 0.5 percent at I < 1000, but by 1 to 10 percent at 1000 < I < 10000. The size of the 
effect exceeds the predicted accuracy of future lensing power spectrum measurements and will likely 
be detected. Precise determinations of cosmological parameters with weak lensing, and studies of 
small-scale fluctuations and clustering, therefore rely on properly including baryonic physics. 

Subject headings: gravitational lensing—dark matter— cosmology: theory— galaxies: formation 


1. INTRODUCTION 

Weak gravitational lensing directly measures the pro¬ 
jected mass distribution and is emerging as one of the 
most powerful and robust probes of the large scale 
structure of the universe, and the nature of dark mat¬ 
ter, dark energy an d gravity llHutereHl2nnl lHll2?inl 
jjjiri fc .1317112004 Ilsba.keTTl] 

2005HKnox et al.ll2005ll . Ongoing and upcoming surveys 
such as CFHTLS\ DES^, LSST^, Pan-STARRS^, SKA'^ 
and SNAP® will significantly reduce statistical errors in 
lensing power spectrum measurements to the <1% level 
at I ~ 1000. As new analysis techniques enable a sig- 
nifica nt reduct i on of systematic errors in cosmic shear 

S is fc .laini 12004 I Jain et ahl 1200,4 iHevinans et al.1 
and references t herein) and cosmic magni¬ 
fication measureme nts llMenard fc Bar telmannI 120021: 
iScranton et~Zll2004 iZhang fc Penll2005alb[l. weak lens¬ 
ing measurement is quickly entering the precision era. To 
match this accuracy, many simplifications i n theoretical 
mode l ing have to be scrutin i zed in detail llBernaidea^ 
19981 iSchneider et alJ 11994 iDodelson fc Zhand 12004 
Dodelson et al.ll200,4 IWhitell200,5ll . 

A widely adopted simplification in weak lensing mod¬ 
eling is the assumption that baryons trace dark matter 
perfectly. With this simplification, weak lensing involves 

^ http://www.cfht.hawaii.edu/Science/CFHLS 
^ http://www.darkenergysurvey.org/ 

® http://www.lsst.org/ 
http://pan-starrs.ifa.hawaii.edu/ 

® http://www.skatelescope.org/ 

® http://snap.lbl.gov/ 


only gravity and collision-less dark matter dynamics, al¬ 
lowing it to be accurately predicted with the aid of N- 
body simulations. However, on small scales, baryons 
do not follow th e dark matter distrib ution. Recently, 
Iwhit^ ll2004l and iZbau fc Kno^ ll2004l estimated the ef¬ 
fects of cooling gas and intra-cluster gas on the lensing 
power spectrum, respectively. These two components of 
baryons can both have an effect of a few percent on the 
lensing power spectrum at I ~ 3000, but with opposite 
signs. Such effects exceed future measurement errors and 
are certainly relevant in order to exploit the full power of 
precision measurements of weak lensing. However, ana¬ 
lytical models are simplified with ad hoc parameters and 
lack the ability to deal with back-reactions of baryons on 
the dark matter. They can not robustly answer several 
key questions such as (1) to what level the two effects 
cancel, (2) at what scales each of them dominates, and 
(3) how large the effect is for the non-virialized inter- 
galactic medium where most of the baryons reside. To 
quantify the effect of baryons on lensing statistics accu¬ 
rately, it is hence necessary to use hydrodynamical sim¬ 
ulations with all relevant gas physics included. In this 
letter, we analyze a controlled set of simulations to ad¬ 
dress this issue. 

2. BARYONIC EFFECT ON THE CLUSTERING OF 
COSMIC MATTER 

2.1. Simulations 

We use the GADGET2 code llSoringel et al.l l2001t 
iSoringell 1200,5jl to simulate structure formation in a 
concordance cosmological model. The cosmological and 
initial density fluctuation parameters of the model are 
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(O™, rJA, ^b, 0-8, n, /i)=(0.268,0.732,0.044,0.85,1,0.71). 
Three simulations were run with the same initial 
fluctuations for a cubic box of 100 h“^Mpc. The 
hrst simulation is a pure dark matter cosmological 
simulation. The second one is a gas simulation where 
no radiative cooling of gas is considered. The last run 
is a gas simulation that includes the physical processes 
of radiative cooling and star formation. It also includes 
supernova feedback, outflows by galactic winds, and 
a sub-resolution multi phase model for the i nterst ellar 
medium as detailed in ISoringel fc Hernauis 3(|2003). In 
all simulations, we use 512^ particles to represent each 
component of dark matter and gas. 

Uncertainties in simulating star formation could bias 
our study. To check the robustness of our star formation 
run, we measured the stellar mass density U* in units of 
the critical density at z = 0. We found U* = 0.0034. 
This is in reason able agreement with r ecent observa¬ 
tional resu lts (e. g. iFukugita et alJ l|1998[ U* = 0.0035), 
IRell et akl ll2nM U* = 0.0028 ± 0.0008)). The star for¬ 
mation history of the simulation is therefore close enough 
to reality for the purposes of the present paper. 

2.2. Clustering in the non-radiative run 

We measure the mass power spectrum P{k) for matter, 
gas and stars separately, as well as for the total matter 
density in the simulations. The Triangular Shaped Cloud 
(TSC) interpolation function W{r) is used to obtain the 
smoothed density held on a uniform grid of 1024^ cells, 
and the power spectrum is determined from the Fo urier 
transfo rm of the density held. We use the method of lJind 
il2?in^ to correct for the smoothing and aliasing effects 
caused by the TSC mass interpolation and for the noise 
caused by the particle discreteness. 

The power spectrum of each matter component is plot¬ 
ted in Figure 1 for redshifts z = 0 and z = 1. Here 
we use the mass variance A^(fc) per logarithmic inter¬ 
val in wavelength which is related to P{k) by A^(fc) = 
4:Trk^P(k)/(2Tr)^ According to the extensive tests of 
iHou et akl (1200511 . the power spectrum is affected by the 
force resolution at A: > w 0.3 x 2Tr/r] , where rj 
is the Plummer-equivalent force softening length used 
in the simulation. In our simulations, r] = 5/i“^kpc 
for the pure dark matter and non-radiative runs, and 
rj = 9/i“^kpc for the star formation run. Down to the 
scale of fc « 20 hMpc”^, the force resolution should have 
negligible effects. 

Comparing the results of the pure dark matter and 
the non-radiative runs, one finds that the gas has a 
weaker clustering than the dark matter on small scales 
(fc > 1 hMpc”^). This can be interpreted as the result of 
gas pressure whi ch reduces small-sca le structure in the 
gas distribution. iZhan fc KnoxI (1200411 approximated the 
hot baryon distribution by assuming that t he gas is in hy¬ 
drost atic equilibrium in NFW dark halos (|Ma,kino et al.l 
[Ml), and found a qualitatively similar result. But quan¬ 
titatively, the gas clustering suppression in their simple 
analytic estimate is stronger than our simulations indi¬ 
cate. At fc = lOfcMpc”^ and z = 0, the A^(fc) of gas 
is 25 percent lower than that of dark matter in our non- 
radiative simulation (see Figure |2Jl, while the A^(fc) of 
gas is 50 percent lower in their analytic estimate. We 
also see a more extended baryonic effect at fc < 1 fcMpc”^ 



k [h Mpc"'] k [h Mpc"'] 


Fig. 1.— The power spectrum of density fluctuations of each mat¬ 
ter component in the non-radiative gas simulation (bottom panels), 
and in the gas cooling and star formation simulation (top panels), 
at redshifts 2 = 1 and 2 = 0. The results are compared with the 
spectrum of the pure dark matter simulation (black solid lines). 
The colored dotted lines, dashed lines, long dashed lines, and solid 
lines are plotted for the dark matter, gas, stellar and total matter 
density fleld, respectively. 


which may be caused by hlaments and is beyond the ex¬ 
ploration of their halo model. These effects indicate the 
need for hydrodynamical simulations to accurately assess 
the effect of baryons on the clustering of cosmic matter. 

An interesting result from our non-radiative simulation 
is that its dark matter has a stronger clustering than 
found in the pure dark matter run. To show the dif¬ 
ference clearly, we plot it separately in Figure [21 where 
we also show the clustering differences of other matter 
components relative to the pure dark matter run. The 
dark matter clustering is a few percent stronger in the 
non-radiative cooling simulation at fc > 1 /iMpc”^, and is 
more than 10% stronger at fc > lOfcMpc”^. This behav¬ 
ior can be understood as a result of the gravitational back 
reaction due to hot gas, since the gas is hotter than the 
dark matter virial temperature in the centr al regions of 
dark halos, due to its collisional nature (see, iRasia et alJ 
l2f)nllbin et'ani2nn^ . 

Compared to the power spectra of gas and dark matter, 
the P^{k) of the total matter distribution in the non- 
radiative run differs less dramatically from the FbM(fc) 
of the pure dark matter run^, because the baryonic sup¬ 
pression of the gas clustering is partially counterbalanced 
by its effect on the dark matter. However, there still ex¬ 
ists a significant difference in the total matter cluster¬ 
ing between the two runs (see Figures 1 and 2). The 
total matter spectrum P^°^{k) is a few percent lower 

^ Here we use the subscripts DM, NR, SF for the pure dark 
matter, non-radiative, and star formation runs, respectively, and 
the superscripts dm, gas, star, and tot for each matter component. 
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Fig. 2. — The influence of baryons on the clustering of each mat¬ 
ter component. The plot shows the relative difference (in percent) 
of the power spectrum of each matter component in the gas simula¬ 
tions relative to the pure dark matter simulation. The results of the 
non-radiative run are plotted with red color, and those of the star 
formation run with green colors. The colored dotted lines, dashed 
lines, and solid lines are plotted for the dark matter and gaseous 
components, and the total matter density field, respectively. The 
stellar component in the star formation run is omitted here, as the 
difference can be easily infered from Figure [Tl 


than Pumik) at k Ki 2/iMpc and is a few percent 
higher than PDM(fc) at fc « lO/iMpc” . At the smaller 
scales, the difference becomes smaller because P^'^ik) 
drops faster than P^^{k) increases. Although the for¬ 
mal gravitational resolution limit of the simulation is 
k w 40hMpc~^, hydrodynamical convergence can only 
be expected on somewhat larger scales, thus it will be 
interesting to check the results at fc > lOhMpc”^ shown 
in Figure 2 with higher resolution simulations. 

2.3. Clustering in the simulation with star formation 

In our star formation run, the stellar matter is found to 
be significantly more clustered than dark matter or gas 
(Figure 1). Its power spectrum Pg^’^ik) is a few times 
higher than the dark matter counterpart Pg™{k), and the 
bias factor increases towards smaller scales. This result is 
expected, because the stars can form only in high density 
regions and thus are highly biased tracers of the matter 
distribution. 

On the other hand, the power spectrum of gas 
is significantly lower than T’g^(fc). Apart from the mild 
effects of shock heating and gas pressure that can be 
infered from the non-radiative simulation, the main rea¬ 
son is that gas is under-represented relative to the dark 
matter in high density regions, where part of the gas 
has cooled and turned into stars. Because of this anti¬ 
bias, P|'p’*(fc) is always smaller by ~ 10% than Pg™(fc) 
even on large scales (Figure 2), and drops much faster at 
small scales than the counterpart P^’^ik) in the non- 


radiative run. The power spectrum of dark matter 
Pg^(fc) looks very similar to the counterpart Pnr(^) 
in the non-radiative run, except the former is higher at 
fc > 10 /iMpc”^ because of the baryon condensation in 
the star formation run. 

The power spectrum of the total matter density Pgp (k) 
shows a similar decrease at fc < 2 /iMpc”^ as found in the 
non-radiative run. As we discussed before, this behavior 
is mainly due to shock heating and the thermal pressure 
of the gas. This is confirmed in our star formation run, 
indicating that this feature is robust against the star for¬ 
mation processes. The Pgp (fc) is higher than P^^{k) at 
fc > 10 /iMpc”^ at z < 1, which is mainly caused by the 
baryon condensation due to gas cooling. We also note 
that PgF*(fc) is slightly lower at fc « 8/iMpc~^ than the 
counterparts P^%ik) and Pouik) at z > 1, which likely 
results from the feedback heating due to active star for¬ 
mation at high redshift. 

3. THE BARYONIC EFFECT ON WEAK LENSING 

In this section, we study the influence of baryons on 
the weak leasing she ar power spect rum Ci. We use Lim¬ 
ber’s approximation llLimberl|1954ll to calculate Ci from 
the simulated pt°*. For a flat universe, Ci and P*°‘ are 
related by 

o = (7) . 

( 1 ) 

where W(x,Xs) = (1 + -z)a(1 “ x/Xs) is the lensing ker¬ 
nel, and Xs and % are the comoving angular diameter 
distances to the source and lens, respectively. 

In Figure 13 we show the effect of the baryons on 
the weak lensing power spectrum by assuming that the 
lensed sources are at redshifts Zg = 0.6, I.O, and 1.5, 
respectively. For the non-radiative run, the lensing 
power signal is suppressed by less than one percent at 
100 < I < 1000 and is then enhanced to about 1% at 
I = 4000. The detailed behavior depends on the source 
redshift, and the relative change of Ci increases with the 
decrease of Zg. This is expected since the baryons have 
more significant effects at lower redshifts (Figure [3). 

Including more realistic star formation processes com¬ 
pensates for the lensing power suppression at 100 < I < 
1000 seen in the non-radiative run. As a result, the 
change of the lensing power by the presence of baryons 
is smaller than 0.5% for Zg > 0.6 and 100 < I < 1000. 
But at 1000 < I < 10000, the baryons can increase Ci 
by up to 10 percent depending on Zg. Again, the relative 
change of C; increases with the decrease of Zg. _ 

Com pared to the previous study of IZhan fc KnoxI 
(I2004^ for the effect of hot baryons, our prediction for 
ACifCi in the non-radiative run is quite different from 
their finding of a monotonic decline of AC; /C; with I at 
I > 1000. We argue that the main reason for this is that 
they neglected the back reaction of the thermal baryons 
on the dark matter in their analytical modeling. 

To match the accuracy of future lensing surveys, mat¬ 
ter power spectrum at scales discu ssed in this Letter 
must be calibrated to ~ 1% accuracy l|Huterer fc Takad'^ 
Thus, the baryonic effect is non-negligible for fu¬ 
ture lensing analysis. Here, we do a simple estimation to 
demonstrate this point. The statistical errors of lensing 
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Fig. 3.— The effect of baryons on the shear power spectrum Q, 
expressed as relative difference (in percent) to the pure dark matter 
model. Here the source galaxies are assumed to be at redshifts Zs = 
0.6, 1.0 and 1.5, respectively. The dashed lines are the predictions 
from our non-radiative run, while the solid lines give the results for 
the star formation run. 


power spectrum measurements (assuming Gaussianity) 
are 


ACi 

Cl 


1 % 


1000 100 0.1 J V ngCi j ’ 


( 2 ) 


where Ug is the galaxy number density, with typical value 
< 100/arcmin^, and 7rms ~ 0.2 is the rms fluctuation 
caused by the galaxy ellipticities. At I < 3000, the shot 
noise term is sub-dominant. For future lensing missions 
with fractional sky coverage /sky ^ 0.1, changes in C/ 
caused by baryons are comparable to, or larger than 
statistical errors, at I > 1000. To constrain cosmol¬ 
ogy using measured Ci’s at all accessible scales (1 less 
than several 10^), the effects due to baryons have to be 
taken into account. Otherwise the derived constraints 
can be biased. For exam ple, since Ci oc (e.g. 

IVan Waerbeke et al.l2001ll . neglecting this effect can lead 
to an overestimate of by several percent. Also, 

it can cause an overestimate of the initial power index 
n. To avoid such biases, an accurate modeling of the 
baryonic effects using hydrodynamical simulations is re¬ 
quired. 

On the other hand, a precision measurement of the 
lensing signal at I of about a few thousand could be used 


to observationally determine the baryonic effects and to 
constrain the galaxy formation process. This can be done 
in two steps. One first uses the lensing power spectrum 
at I < 1000 to constrain 12^, f^DE, crs, etc. Even though 
this uses only part of the information in the lensing power 
spectrum, the cosmolo gical constraints obtaine d are not 
significantly degraded liHuterer fc Whit jl200fill . The de¬ 
rived cosmological model is then used to predict Ci (and 
its statistical fluctuations) at I > 1000, assuming no 
baryonic effects. One can then quantify the deviation 
of the measured C/ at Z > 1000 from the predicted C/. 


4. CONCLUSIONS 

In this Letter, we used a controlled set of high resolu¬ 
tion N-body and hydrodynamical/N-body simulations to 
study the influence of baryons on the clustering of cos¬ 
mic matter. Since the three simulations we used have 
identical initial condition, simulated power spectra suf¬ 
fer essentially the same sample variance. Since we quan¬ 
tify the baryonic effect as ratios of corresponding power 
spectra, the result presented in this Letter is effectively 
free of sample variance. In both the non-radiative sim¬ 
ulation and the simulation with gas cooling and star 
formation, the clustering of the gas is suppressed while 
that of dark matter is boosted at fc > l/iMpc”^. The 
stellar mass is highly biased relative to the dark mat¬ 
ter in the pure N-body simulation. Despite a partial 
counterbalance between the dark matter and the gas, 
the clustering of the total matter is suppressed by up 
to 1% at 1 < fc < 10hMpc“^, and is boosted up to 
2% in the non-radiative run, and 10% in the star for¬ 
mation run at fc w 20/iMpc“^. Using these power spec¬ 
trum measurements to study the baryonic effect on the 
weak lensing shear measurement at 100 < I < 10000, 
we found that baryons can change the shear power spec¬ 
trum by less than 0.5% at I < 1000, but by 1% to 10% 
at 1000 < I < 10000. Therefore, the influence of baryons 
on the clustering of cosmic matter will be detected in 
future weak lensing surveys. Understanding these bary¬ 
onic effects is not only important for galaxy formation, 
but also crucial for accurately determining cosmological 
parameters with cosmic shear, and for constraining the 
initial fluctuations on small scales. 


The simulations were run at Shanghai Supercomputer 
Center. The work at Shanghai is supported by the 
grants from NSFC (Nos. 10125314, 10373012, 10533030) 
and from Shanghai Key Projects in Basic research (No. 
04JC14079 and 05XD14019). PJZ is supported by the 
One Hundred Talents project of Chinese academy of sci¬ 
ence. 
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